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A complete, non-demolition procedure is established for 
measuring multi-qubit entangled states, such as the Bell- 
states and the GHZ-states, which is essential in certain pro- 
cesses of quantum communication, computation, and telepor- 
tation. No interaction between the individual parts of the 
entangled system, nor with any environment is required. A 
small probe (e.g. a single qubit) takes care of all interaction 
with the system, and is used repeatedly. The probe-qubit in- 
teraction is of the simplest form, and only this one type of 
interaction is required to perform a complete measurement. 
The process may be divided into elementary local operations 
and interactions, taking place sequentially as the probe visits 
each of the qubits. A shuttle mode is described, which may be 
repeated indefinitely. By the quantum Zeno effect, the entan- 
gled states can be maintained until released in a predictable 
state. This shuttle process is stable, and self-correcting, by 
virtue of the standard measurements performed repeatedly on 
the probe. 

PACS Numbers: 03.65.Bz, 89.70,+c, 42.50.Dv 



I. INTRODUCTION 

Quantum entanglement is one of the most remarkable 
manifestations of the fundamental principles. Consider 
two distinguishable systems, with quantum states | ipi ) 
and \(f>i), i = 1, 2 • • •. When these systems are indepen- 
dent, then their joint state is separable, as for instance 
I tpi ) <8> | (f>i ) , or | ?/>2 ) ® | 02 )• However, according to 
the principle of superposition, the compound system may 
also exist in states that are not separable, such as 

I * > = Ci| lp! ) <8> | 01 ) + C 2 | V>2 ) ® I 02 > + ■ ■ ' 

where cx, C2,- • • are complex amplitudes. This possibility 
has given rise to many profound discussions, in particu- 
lar concerning the EPR paradox and the Bell inequali- 
ties . It has been confirmed — beyond any reasonable 
doubt — that these entangled states are uniquely a quan- 
tum phenomenon |j| |5| . 

There have been many proposals for preparing specific 
entangled states, starting with suitably initialized states, 
such as for instance [pMllI] , and this has been achiev ed in 
new ways, for Bell-states |f 2|-|l5f (notation, cf. Sect. [I A ) 



h(\ TT)±IU» 



V2 



\^) = ^(\ U>±IIT» 



and for GHZ-states 10-2 



|s ± ) = ^(|TTT)±IIU» 

Recent experiments have demonstrated that it is feasi- 
ble to design processes that depend in an essential way on 
handling such entangled states. This includes quantum 
communication pi) , quantum computation ]22|-p8||, and 
quantum teleportation |2^-|3^]. For instance, in densely 
coded quantum communication the receiving party must 
measure the BeU-states in order to retrieve the encoded 
information p^ ]. In quantum teleportation [ p5[ the ex- 
ecutive step is also a measurement of the entangled BeU- 
states, or an EPR-state p6| , p7[ . This projects the en- 
tire system into the teleporting configuration, and, at the 
same time, acquires the data that must be transmitted 
by classical means. 

In order to perform such a measurement one has to 
design a global experimental situation for the compound 
system. For pairs of photons, interferometric methods 
can provide characteristic detection patterns, which more 
or less completely distinguish the entangled states from 
each other |2^,^l|,^8 - 42 1 . Unless interaction is possible a 
complete Bell-state measurement is not feasible p3| , ^4| . 
In |3(| , where the entangled components are the polariza- 
tion and momentum of a single photon, all the BeU-states 
can be distinguished by using polarizing beam-splitters 
|Q. In this situation the entangled subsystems (i.e. po- 
larization and momentum of a single photon) are made 
to interact directly with each other. That is also the case 
in a recent NMR based teleportation experiment |Q . In 
these procedures, what amounts to a controlled-NOT op- 
eration disentangles the Bell-states jl6 47 1, and one can 
then measure the now separate subsystems. The entan- 
gled state is destroyed in these types of measurement. 

However, states of known entanglement are becoming 
a resource. It will be desirable to retain the entangled 
system itself after the measurement, for subsequent pro- 
cessing — not merely the data acquired. 

Also, from a more fundamental point of view, in order 
to establish entanglement as a standard physical observ- 
able, one really needs a non-demolition procedure. That 
is, a procedure which leaves any of the properly entangled 
states invariant, while delivering unambiguous informa- 
tion about it. A method of this nature has been proposed 
for the BeU-states and other entangled states of a similar 
structure in |36],[l8) . In this design, the measuring device 
itself consists of components that must be prepared be- 
forehand into entangled states. After these components 
have been placed at the relevant locations, an instanta- 
neous non-local measurement can be performed, at least 
in principle. 



1 



The measurement procedure established in the present 
paper is different from the existing ones in several re- 
spects. It is both simple and economical of resources 
that appear to remain scarce, at least within a fore- 
seeable future. The 'apparatus' can be as small as a 
single-qubit probe, which is used repeatedly according to 
a predesigned shuttle schedule. In principle, this allows 
a complete measurement of the entangled basis-states of 
any number of target qubits. The necessary interactions 
have been reduced to bare essentials, so that only the 
most elementary operations are involved at every stage 
of the procedure. In particular, no entanglement is re- 
quired at the outset. For example, it is possible to imag- 
ine a single 13 C nuclear spin probe, paying simultaneous 
attention to a number of different H qubits. 

Since the method relies on a traveling probe, it con- 
forms in a straightforward way with special relativity. 
On the other hand, this means that it is not capable 
of instantaneous measurement at space-like separation. 
However, it seems that it could be efficient in combi- 
nation with entanglement swapping, in creating distant 
multi-qubit entangled states 

Sect. H describes the principles of the procedure for 
measuring the Bell-states of two qubits. Examples in- 
clude both a spin 1/2, and a single- mode cavity-field 



probe. Sect. Ill 



entends these methods to systems of 
arbitrary numbers of qubits, and provides an extensive 
analysis of the GHZ-states of three qubits. For both cases 
a shuttle design is provided, which permits the measure- 
ment to be continued indefinitely, with a sustained flow 
of measurement data. In effect, this could provide stable 
storage for entangled states, which can then (in principle, 
of course) be extracted with certainty at a predetermined 
time, using the acquired data. 



II. MEASURING THE BELL-ASPECT 

A. Definitions 

Consider two qubits, labeled 1 and 2, and represented 
as distinguishable spin 1/2 systems with Pauli spin op- 
erators a = (a x ,a y , a z ). Let a z \ | ) = | T > and a z \ J. ) = 
— 1 1 ), and write their joint states in the form 



U)^IT) (1) ®II) 



(2) 



The Bell-states form an orthonormal basis of (maximally) 
entangled two-qubit states 



|* ± ) = ^(|U>±UT)) 



|$ ± ) = -i 5 (|H)±IU)) 



(2.1a) 
(2.1b) 



These states are simultaneous eigenstates of the commut- 
ing operators 



XX — u x 

yy — u y 

.... = rrW 



r (2) 

X 

r (2) 
V 

J2) 



(2.2) 



with eigenvalues ±1. They are therefore also eigenstates 
of the 'Bell-operators' ffiM 



-BCHSH 



2(<j xx ± a 



mil 



The eigenvalues ±2V2 signal a maximal violation of the 
corresponding Bell- inequalities, in the well-known way. 
The combination of this set of compatible observables, 
and their basis states, will be referred to as the 'Bell- 
aspect' ]50[ . 

The eigenvalue of a xx determines the 'e/o question', 
i.e. whether the state is even or odd (e/o) under global 
spin flip. The eigenvalue of a zz determines the 'p/a ques- 
tion', whether the spins are parallel or antiparallel (p/a). 
The corresponding projectors are 



Ph 



a = ±1 
b = ±1 



(2.3) 
(2.4) 



Consequently, the Bell-aspect is defined by the one- 
dimensional projectors 

(2.5) 



Pab = |(l + aa xx )(l + ba zz ) 



A complete measurement can therefore be done in two 
stages. Each stage consists of a partial measurement sep- 
arately deciding the p/a and e/o questions. Th e rel evant 
procedures must commute, as the projectors fl2.3| ) and 
(0 do. 

Let a bilateral rotation of the two spins, by —tt/2 about 
the global y-axis, be written as 



Uyy Uy ® Uy , U y 



V2 



1 + ier,, 



Then 



(Jx = U\a z U yi (Uyf = i(T y 
Using the algebra of the Kronecker product one gets 



— ^yy° zz ^yy 

= -(Uyyf = 



Therefore 

Pah = 



-jcr y yUyy(l + aa zz )U yy (l + b(7 zz ) (2.6) 

This suggests a simplified two-stage procedure, consisting 
of identical operations U vy (1 + aa zz ). Here one first de- 
cides the p/a question, then one rotates the spins bilater- 
ally. With such a chain of operations, the measurement is 
reduced to the most economical form, with respect to the 
resources that will be needed to carry it out in practice. 
In particular, only one type of unilateral rotation, U v , 
is required. The essential task is to perform the partial 
p/a measurement, but also here only a single procedure 
is necessary. What is more, each partial measurement is 
a binary test. It therefore requires no more than a single 
spin 1/2, repeatedly probing the qubits. 
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B. Spin 1/2 probe 

Let the probe eigenstates in the xy-plane be denoted 
IE), |W), IN), and |S), i.e. 



crJE) 



IE) 

IN) 

ct 2 |E) =|W) 



a v \N) 



o-jw; 



|W) 



a z |W) = |E) 
The corresponding projectors are written 



P E 

P N : 



|E)(E| 
N)(N| 



Pw 
P S = 



|W)(W| 
S)(S| 



(2.7) 



In the following, a — (cr x ,a y ,a z ) (with no superscript) 
stands for the Pauli operators of this probe. The explicit 
<g> separates the probe and the qubit operators. 

The interaction between a single qubit and the probe 
is taken to be of the simplest possible form 

Ux=exv{-i9a z ®a£ ) /2) 

The parameter can be adjusted by controlling the 
length of time during which the probe interacts with 
qubit 1. This type of interaction can be realized in differ- 
ent systems, including nuclear spins, where it is standard, 
and by dispersive Rydberg-atom/ca vity-fi eld interaction 
(described in the following, cf. Sect. II H ). It allows the 
entire measurement process to be reduced to elementary 
operations. This strategy is of course familiar from NMR 
work, for instance, which establishes its feasibility. 
A similar interaction with qubit 2 gives 

U 2 =cxp{-i9cr z <g>a^/2) 

It does not matter if these interactions with the qubits 
take place simultaneously, or one after the other. All 
operators commute, and 



U 2 Ui = exp(-i6»cr 2 (g) S z ) 



where 



(2)\ 



Extend the definitions (2.2) to include the unit operators, 
such that for instance 



1 



(2) 



Then one can use aoz — Czo^zz to write 

S z = \ (czo + foz) = cr z0 P+ = CT0Z-P4 



(2.8) 



The projector is given by (2.4). For j = 1,2, 



VzoP+f 3 

(<J Z ® <JzoP+f J+1 = (a 
Therefore 



10P-4- 



icr z o)(l®P + 



U 2 U X = 1 <g> 1 + (cos0 - 1) (1 <g> P+) 
-ism9 (a z <g> er z0 ) (1 <8> P+) 
= l(g>P- (2.9) 
+ (cos 9 — ?sin# a z ® a z o) (1 ® P+) 



For 9 = vr/2 then 

U 2 Ux = 1 ® P- - ia z ® <j z0 P+ 



(2.10) 



Here a z rotates the probe spin, while <j z q swaps the qubit 
Bell-states, e.g. 

a z0 \^) = \^) 

Because of the P+ projector, both of these operations 
take place only when the qubit spins are in the parallel 
subspace spanned by $ . 



C. The partial measurement (p/a) 

Let the initial states of the probe and qubits be given 
by statistical operators (density matrices) v and w, re- 
spectively. Suppose the probe is initialized along the pos- 
itive x-direction, so that the joint state is 

p = v ®w, v = Pe 

After the probe-qubit interactions one has 

P ' = u 2 u lP u\ul 

= |E)(E|(g)P_wP_ 

+ |W)(W|®a z0 P + u>P+(7 2 o 
+i|E)(W|® P-wP+cr z0 
-i I W)(E I ®a z0 P + wP- 

One can now read the probe by any ordinary, local mea- 
surement, coupling it to a suitable apparatus, possibly 
including a large environment. As far as the probe and 
qubit degrees of freedom are concerned, this is equivalent 
to the operation 



Pmcas = (Pe ® 1) P (Pe <8> 1) 

+ (P w ® 1) p (Pw <g> 1) 



(2.11) 



That is, none of this should affect the qubit pair. Mea- 
suring the probe spin in the x-direction then gives 

Pmcas = Pe <8> P-wP- + P w <8> <J z0 P + wP + cr zQ (2.12) 

This represents the acquisition of binary data, 'E' or 'W 
corresponding to the distinct alternatives of the p/a test. 
The (reduced) state of the qubit-pair is not affected by 
these interactions between the probe and the external 
apparatus 

w' = tr pro be(p / ) = tr pro bo(Pmcas) 

= P_wP- + a z0 P + wP + <7 zQ 
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The statistics of the data is that, one obtains the item E 
a fraction of the times given by 

tr((P E ® l)p'(P E ® 1)) = tr 12 (P_wP_) 

= tri2 (wP_) = p- 

This equals the estimated probability, p_, for anti- 
parallel spins in the initial qubit state w. Likewise, one 
gets W with a frequency 

tri 2 (cr z oP+wP + cr 2 o) = tri 2 {wP + ) = p + 



D. Probe motion 

It is also of interest to monitor the probe spin in the 
course of time during the interaction with the qubits. 
The reduced probe state at a time corresponding to 9 is 
given by (using (|2.9|)) 

v(6)=tr 12 (U 2 U lP uM) 

= (p- + p+ cos 2 0) v + p + sin 2 9 o z va z 
+i sin 9 cos 9 (S z ) (va z — a z v) 



Here, according to (2.8) 



(S z ) = tr 12 (wa z0 P + ) 

= (TT|u,|TT>- 



($+ | w |$" 

;umu) 



This quantity is zero in all the Bell-states, but not in the 
triplet states where S z = ±1. The probe spin polariza- 
tion is 

m(0) = trprobc^fl) a) 

This gives 

m x {9) = (p- +p+ cos26») m x (0) - sin 29 (S z ) m y (0) 



m v {9) = (p_ +p+< 
m z (9) = m z (0) 



26>) m y (0) + sin 2(9 (£*) m x (0) 



For a large ensemble this motion may be detectable, 
and would then provide information about w in terms 
of the mean values p± and (S z ). If no probe measure- 
ment is made, then (disregarding decoherence) every- 
thing eventually returns to the initial state p. On the 
other hand, in the measurement situation, where one 
starts at rh = (1, 0, 0), one gets instead at 9 = tt/2 



■ P-Pe+p+Pw 



and 



m x (-K/2) =p--p + 
m y (-K/2) = to 2 (tt/2) = 

Of course, a measurement on the probe at this stage re- 
moves its coherence with the qubit pair, as it should, 



TABLE I. Bell-state propagation during shuttle operation 
(phases omitted). The probe starts in the state Pb- 
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and the probe-qubit system therefore cannot return to 
the initial state. The data thus acquired determines the 
initial state of each individual probe-qubit system in the 
subsequent evolution. In the case of W this means that 
the probe will be in the state Pw, and, at the same time, 
the qubit pair state is projected into the parallel sub- 
space. This occurs with a statistical frequency equal to 
the probability p + . 



E. Repeated measurement 

One can cancel the $ ± permutation due to a z o by per- 
forming a repeated measurement, using the output as 
new input. Then, using (2.12) as initial state, either 



P E <g) P-wP- h-> P E ® P_wP- + 

or 

Pw ® <J z0 P + wP + a z() 

h-> + P E <£> a z0 P+(T z oP + wP + a z0 P + a z0 
= P E (8 P+wP+ 

Assuming one starts the probe initially in P E , this pro- 
duces data EE and WE, respectively, and the probe al- 
ways ends up in P E . The qubit pair experiences a clean 
projection 



Repeated — Pe ® (P-twP- + P + ^P+) 



(2.13) 



Although somewhat extravagant, this duplication of the 
data is of course a useful check, and the resulting state is 
the ideal input for the second stage in a complete mea- 
surement. It is essential that one perform the probe mea- 
surement at the end of both of the U2U1 interaction se- 
quences. Another way to compensate for the unitary 
transformation by er z o will be described in Sect. Ill C| . 



F. Second partial measurement (e/o) 



According to (2J3) a complete measurement consists of 
a sequence (reading from right to left and omitting the 
final a yy ) 



{M x U yy U 2 U 1 ){M x U vy U 2 U l ) 



(2.14) 



Here, M x stands for any (unspecified) interaction of the 
probe with its supporting apparatus, so that the outcome 
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is the required probe measurement in the x-direction. 
There are many different ways to do this, and it is only 
nec essar y that the result is the measurement operation 
of ( 2.11 ), as far as the probe's state is concerned. The 
simplest presumably is a process like the one analyzed in 
|plf . Such a dephasing has recently been demonstrated 
in NMR teleportation [[33| , taking advantage of an order- 
of-magnitude difference in the T\ and T2 time-scales on 
the 13 C spins. 

The first measurement decides the p/a test, the sec- 
ond one the e/o test. The bilateral rotation U yy of the 
qubits commutes with any M x . As mentioned, since the 
probe is projected into either | E ) or | W ) , it is automati- 
cally initialized for the following measurement. Aft er th e 
second measurement the state has become (using ( j2.12| )) 



p" = Pp 

rmcas & 



{UyyP-UyyP-)w(dO. 



+P W ® (U yy a z0 P + U yy P-) w (do.) f 
+P W <E> {U vv P_Uy V a z0 P + ) w (do.) f 
+Pe ® (U yy a z0 P + U yy a z0 P + ) w (do.) f 

Here U yy causes the following permutation 

tf W |¥+) = |3->, Uyy\ 

and leaves the other two Bell-states invariant. The 
handling of the Bell-states is summarized in Table Q. 
ft is understood that, permutations between the basis 
states within the Bell-aspect are acceptable (i.e. non- 
demolition). Else they can be compensated. 

More generally, for any input state, w, of the qubit sys- 
tem, one is now in possession of pure Bell-states, which 
emerge with statistical frequencies equal to the theoreti- 
cal probabilities predicted with that w, i.e. 

Pab = tTl 2 {Pabw) 

To verify this, consider that any mixed state can be writ- 
ten as a convex combination of pure states 




(b) 




FIG. 1. Shuttle-operation in the complete, non-demolition 
measurement of the Bell-aspect. Equivalent diagrams, for a 
traveling probe (a), and a sta tiona ry probe (b). The wiggly 
lines represent the interaction (2.10) between the probe p and 
the qubits. U y are local —90° qubit rotations, and M x are 
standard spin 1/2 measurements, on the probe exclusively. 

of cases, equal to p- = |c + _| 2 + |c | 2 , where the probe 

is not flipped, this produces 

|Ve> = (c+_|$-) + c__|*-))/VpT 

Otherwise, the probe is flipped, and one gets 

I ihr) = (-c ++ | *+ > + c-+| <&+ » /VPT 

The second measurement now gives | ipEE ) = | ) hi 

a fraction of the first cases equal to |c \ 2 /p-- Likewise 

for the alternative cases. This confirms that one obtains 
the pure output states listed in Table Q with statistical 
frequencies equal to 



The are weights (i.e. Yl £7 = l>£j — 0)> an d the pure 
states I ipj ) need not be orthogonal. Any pure state of 
the qubit pair can be expanded in the Bell-basis 



|V)=c++|$ + ) 
where 



E i c ^i 2 = 1 

a,b=± 

During the first measurement, there is a projection of | ip ) 
into the subspaces of parallel and anti-parallel spins, and 
a subsequent bilateral rotation by U yy . In that fraction 
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Pab 



G. Shuttle mode 

All interactions are local. It is therefore possible to op- 
erate the present procedure in shuttle mode, where the 
probe can travel between qubits located in different re- 
gions of space. Of course, an instantaneous measurement 
at space-like separation is not feasible by this method. If 
the rotations U y can be done at each qubit site, and if 
there is an M x station, then (2.14) can be arranged as 
follows 
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{UyU 2 ) (M x UyU 2 ) ■ ■ ■ {UyU^ (hlitAk) 



qubit 2 qubit 1 

• • • (UyUx) {M x UyU!) • ■ ■ 

qubit 1 

This is represented graphically in Fig. |l|. 

In this shuttle procedure, both of the probe states Pe 
and Pw occur as input, and what is essential in the data 
is whether the probe was flipped or not. Let '1' stand for 
a flip, and '0' for no flip. There are two types of evolution, 
provided the shuttle action continues. Either, the state 
falls into | ) and there are no flips of the probe. Or, 
there is a cyclic shift among the remaining Bell-states 



After the first roundtrip of the probe, provided its initial 
state is known, the measurement is complete. The data 
uniquely specifies the measurement result with respect to 
the input state. This is what is primarily of interest in 
quantum communication and teleportation. 

One also knows the specific Bell-state at the output. 
Each output Bell-state is identified by the three most 
recent data items, i.e. in terms of probe flips 
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In this cycle, after the second roundtrip one has had 
the opportunity to extract any of the three alternatives. 
That is, after the two first probe measurements one has 
a pure Bell-state, and subsequently every measurement 
presents the next one. Therefore, the shuttle procedure 
may be an economical way of preparing specific Bell- 
states. The only statistical element remaining in this 
preparation procedure, even with completely arbitrary 
input states, for both probe and qubits, is whether the 
system falls into one or the other of the two cycles. 

Another interesting possibility is to employ the stabil- 
ity of the present shuttle cycle for storage. Any perturba- 
tion, or imperfection, will be corrected after one (error- 
free) roundtrip (2 bits of data) . With good quality in the 
operations Ui%, U y , and M x , it would be expected that 
imperfections only result in recognizable glitches in oth- 
erwise perfect cycles. This would allow persistent Zeno 
monitoring of the Bell-state of the qubit pair. If the 
probe-qubit interaction is on continuously, what is re- 
quired is a regularly timed sequence of 'instantaneous' 
M x U y operations. Such a time-scale gap, apparently, 
seems to agree with the given conditions in some of the 
relevant systems. 



H. Cavity-field probe 

Suppose a cavity mode interacts dispersively for a spe- 
cific length of time with a passing Rydberg atom p2[ . 
The resulting transformation is of the form 

This interaction is suitable for the present purpose, and 
has the advantage that it conserves the photon number. 
The number operator n generates a phase-shift in the 
field that occupies the cavity. Thus for two qubit-atoms 

U2U1 = exp(~i26h® S z ) 

= 1 <e> 1 

+ (cos(20n) - l) ® P + - i sin(26>n) ® a z0 P + 
= 1®P- 

+ (cos(2#n) ® 1 - i s'm(26h) ® a z0 ) (1 <E> P + ) 

First, consider using cavity number states | 0) and 1 1) 
to create phase-states ||| 

|E) = ^(|0) + |1)), |W> = ^(|0>-|1» (2.15) 

To achieve phase-flips between these states the interac- 
tion time must be adjusted so that 8 = ir/2. This gives 



U 2 Ui = I® P- + (-l)"(g)P+ 



(2.16) 



One may then consider passing the qubit-atoms through 
two such probe cavities in succession, with a U y rotation 
in between. 

According to the Jaynes-Cummings model at e xact res- 
onance J54|, the cavity field phase-states ( [2.15] ) can be 
transferred to an auxiliary probe-atom, whose state is 
then measured |5^J5f|. The cavity can be initialized in 
an analogous way fl57j], similar to the one used in Jlifl to 
prepare entangled atom states. Another way to measure 
the phase in two opposite directions has been proposed 
in (58). 

The data obtained by measuring the cavity-field phases 
after both atoms have passed would be analogous to those 
obtained with the spin 1/2 probe. With two distinct 
probe systems, it is possible to postpone acquiring the 
data, since there is no physical interaction with the qubits 
involved in this. The qubit state is the mixed 



w„ 



E 

,,b=±l 



PabwPab 



One needs the probe data to separate the four parts in 
this mixture. There is no permutation of the $ states, 
so by including the final U£ y one gets a clean projection 
on the Bell-aspect. If it is omitted, then there is a per- 
mutation corresponding to U yy . 

Next, consider starting with a coherent state of the 
probe field 



.ATT 



n=0 



G 



where | n ) are the number eigenstates, and a a complex 
number. For a qubit state | ip ) one gets 



U 2 Ux\a)®\il)) = \a)®P-\i/>) 



Although the opposite-phase coherent states are not or- 
thogonal, their overlap decreases with field amplitude 



These field states resemble classical, i.e. distinct, pointer 
positions. Furthermore, it has been demonstrated explic- 
itly that the coherences in the field density matrix d ecay 
in time, more rapidly for large amplitudes \oc\ [pqjsgf . 
The time dependence appears to agree with the model 
pPf where the cavity field is coupled to a reservoir, which 
causes amplitude damping: aexp(— At/2), where A is a 
coupling constant, and t the time. According to this 
model, the ratio of off-diagonal elements to diagonal ones 
decreases with time as 



M - a , 

(a I a) 



-exp(— Xt) 



= e 



-2|a| 2 At ; for Xt-^l 



A dispersive coupling allows the qubit-atoms to inter- 
act simultaneously with the field probe. For moderately 
large amplitudes this could, in principle at least, take 
place before any significant decoherence has occurred. 
Subsequently, the probe measurement would be approx- 
imated by the decay of the coherences, still leaving time 
to read the data, before the amplitude reaches zero. Re- 
cently a model has been proposed |3l| in which the field 
amplitude decay is prevented, although there is diffusive 
broadening. 



III. MEASURING GENERAL ENTANGLED 
ASPECTS 

A. The GHZ-aspect 

Of course, eigenstates of a zz , such as | |T ), or a xx , such 
as | —»—►), are not necessarily entangled. But the simul- 
taneous eigenstates, the Bell-states, are. For instance 



•) = IU> + UT) 



(3.1) 



where 



and 



are eigenstates of a x . A three-qubit 



generalization of (2.5) therefore is 

r X xx)(l + ba zz0 )(l + ca zz) (3.2) 



Pabc 



i 



with a,b,c — ±1. The partial projectors commute, since 
for example 

CxxxCzzO — CxxCzz ® Cr x l = O zz O xx ® Id : x = U zz 0Cr xxx 

The 8 entangled states corresponding to P a b (cf. Table |J) 
are mutually orthogonal, since at least one of the quan- 
tum numbers a, b, c is different between any pair. The set 



TABLE II. Basis states of the GHZ-aspect for three qubits. 
The two representations are equivalent, the second alter- 
native using Bell-states of qubits 1 and 2, together with 
o-^-eigenstates of qubit 3. 



abc 



standard basis 



w< 12 >®w< 3 > 



1 


+ + + 


1 TTT) 


+ 1 TTT) 


|$+^ 


+ 


<-) 


2 


+ + - 


TTI) 


+ 1 TTT) 


$+ -» 




$- <-> 


3 


H h 


1 Til) 


+ 1 TTI ) 


*+ -» 




*-) 


4 


H 


TIT) 


+ 1 TTT) 


*+ -> 


+ 


«-) 


5 


- + + 


1 TTT) 


- 1 TTT) 


1*- — 


+ 


$+ <-) 


6 


- + - 


TTT) 


- 1 TTT) 


[*--> 




$+ <-) 


7 


h 


1 TTT) 


-1 TTT) 


*- -> 




*+ ^> 


8 




TTT) 


- 1 TTT) 




+ 


*+ ^> 



will be referred to as the 'GHZ-aspect', since it is built 
to contain the original GHZ-states |Tc] 



l 



(I TTT} ± I III}) 



Obviously, the GHZ-basis states can be turned into each 
other by means of one or more unilateral spin flips, gen- 
erated by a x . Measurement of this aspect was suggested 
for the 9-qubit error correction code |6^| . These states 
also remain of great interest with respect to the purpose 
for which they were originally proposed. 

Again, since one starts with just pairwise p/a tests, 
entanglement requires the global e/o test of the projector 



(3.3) 



The complete, non-demolition measurement proc edure 
for the GHZ-aspect will be presented in Sect. HID. The 
following is a general result, providing modules that per- 
form partial measurements for this and other aspects. 



B. The partial measurement module 

The partial measurement corresponding to the generic 
projector 



+ aa z ... z ) 



is an operational component in a complete measure- 
ment of an n-qubit entangled aspect. It is convenient 
to label the k qubits involved (k < n) sequentially as 
j = 1,2, ... ,k. The remaining n — k qubits are ignored 
during the present procedure. For each individual qubit 
the quantization axis can be determined by local rota- 
tions, combined into a multilateral operation, such as 



This makes the measurement suitable for testing 



any fc-qubit spin correlation. In particular, 

|(l + aa x .. 



Pa — Uy...yPaU>i 



V---V 



The complete measurement of a GHZ-class entangled as- 
pect can therefore be carried out using (essentially) only 
the present type of partial measurement operation. 
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Consider a suitable probe property, A, such as for in- 
stance o~ z /2 or n. Suppose the interaction of the probe 
with the j'th qubit is 



Uj = e - i0m °*' = U(9) $ + U{-8) <g> P w (3.4) 



where 



U{9) = 



iOA 



(3.5) 



As usual, 9 is an adjustable angle, determined by the 
duration of the interaction. It is chosen (since that is 
adequate) to be the same for each j. In order to achieve 
a measurement it is also necessary to choose A, and the 
probe state, carefully. The probe can interact either si- 
multaneously, or sequentially, with the k qubits. 

a. Theorem After these interactions the joint probe- 
qubit system has been transformed by the following uni- 
tary operator 

U k ---U x (3.6) 

= U(a k ) ® ( 7 j 1 - fc) Pj 1 - fc) + U(f3 k ) ® U?- k) P< l - k) 

where 



a k = fcf, (3 k = a k - 7r 



(3.7) 



The opera tor to the left of the explicit ® is for the probe, 
given by ( |3.5[ ). Those to the right are for the k qubits. 
The projectors are (J = 1, • • • , k) 



with notation 



= 1®. 

= (i) 



• ® 1 



The operators L/± are given by 



jj(i-i) = jj(i-i-i)p(i) + ^(i-j-i)pO') 



(3.8a) 
(3.8b) 



6. Proof The proof is by induction. As stated in (^J), 
for the first qubit-probe interaction, 17%, the expression 
agrees in form with ( |3.6[ ). Given that (3.6) is true for the 
first j qubits, then 



[u{9) ® p} j+1) + C/(-6i) ® P 0+1) } 
< {U(aj) 



(3.9) 



(1-3) p(l-3) 



tftft) ® cr 



(l-3)p( 



;l-3) j 



(qubit unit operators are not displayed). One needs the 
identity 



p(i+i) p(l-3) _ p(j+i) p(l-J'+i) 



(3.10) 



Here, on the left-hand side of ( 3.10 ) 

(i + a ^+ 1 ))(i + 64 1 .7 ) ) 



= 1 



r (3+l) 



(3+1) 



where a unit operator has been inserted, an d th is agrees 
with the right-hand side. The four terms in ( |3.9[ ) are then 
given by 

(1-3) D (3'+l)p(l-3+l) 



U(otj + $) ® ift ' Jj P 
t/(ft + 8) ® ^-^pW+^pl 1 -^ 1 ) 

- 0) ® cr(i-i)pCj+i)p(i-i+i) 

f/(^. _ ^(i-i)p(j+i) P _(i-i+i) 



If one takes 9 = n/2, and if ( |3 . 7|) holds up to j, then the 
parameters in the probe U are related as 

«i + 1 = (i + !)! = 

ft + f = (i + 1) f - 7T = ft +1 
«3 - f = C? + 1) § " 7T = ft+1 

ft - § = (.? + 1) f - 2tt = a j+1 - 2vr 
The last can give rise to a phase factor 



r\ = U(-2<k) 



AI-nA 



The recursion starts with 

= U W = 1« 

The projectors in these expressions are for a single qubit 
(the j'th), and 77 is a phase factor to be defined in the 
following. Note that, if 77 = 1 this immediately implies 
that (all j = 1, • • • , k) 

jj^-j) = = 1(1-3) 



The measurement design requires this to be a number for 
all relevant probe states. Evidently, in order to obtain a 
measurement transformation, there has to be restrictions 
of some kind, else one would have just any transforma- 
tion. For A = a z /2 one has 77 = — 1, and for A = h 
one has 77 — 1. If this condition (that 77 be a number) is 
satisfied, then the proof is complete. 

c. Check With a spin 1/2 probe and k — 2 one gets 



a 2 



ft 







and 



u{ X2) = u{ l) p (2) u{ _X )p (2) = p{ 2) p (2) = a ( 2) 



jj(12) = p(2) 



,pf = 1 
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So 



U 2 U 1 = U(tt) ® p} 12) + C/(0) ® f! 12) 



in agreement with (2.10). 

Since the probe advances, or retreats, by tt/2 in each 
encounter with a qubit, all the four compass states (2.7) 
can occur. One gets 

p' = U k ■ ■ ■ U t v <g> w U\ ■ ■ ■ Ul 

= U(a k )vU(a k )^ui h -- k) Pi 1 --- k) wPl 1 --- k) ui 1 ---^ 



TABLE III. GHZ-state propagation during the uncompen- 
sated measurement according to (3.2) (overall phase omitted). 
Theprobe starts in | E ). Numbers refe r to the listing in Ta- 
ble O, letters to the four compass states (2.7) of a qubit probe. 



input 


U 2 Ui 


M x 


U2U3 


M x 




U3U2U1 




1 


5 


W 


1 


E 


|*+ 


-f)-|$- <-) 


s 


2 


6 


W 


6 


W 


|$- 


-►>- |*+ <-) 


s 


3 


3 


E 


7 


W 


j*- 


->) - |$+ <_) 


s 


4 


4 


E 


4 


E 


[*+ 


->}-!#-<-) 


s 


5 


1 


W 


5 


E 


$~ 


+ «-) 


N 


6 


2 


W 


2 


W 


[¥+ 


->) + 1$- <-) 


N 


7 


7 


E 


3 


W 


|$+ 


-►) + !*- <-) 


N 


8 


8 


E 


8 


E 




->) + 


N 



With even k one must measure the probe in the x- 
direction, using M x . For odd k in the y-direction. Of 
course, one can then consider rotating the probe by 7r/2, 
in order to be able to use the M x equipment repeatedly. 

For example, with even k, starting with v — P&, after 
the probe has been measured, off-diagonal terms in the 
probe part of p' meas are °f the form 

PE,wU{a k )P^U(a k - n)^P E ,w 

They will vanish due to the extra rotation by 7T, pro- 
vided the relevant probe states are orthogonal (by de- 
sign). Likewise, the diagonal terms are 



PE,wU(a k )P E U(a k yP E , 



and 



Pw,EU(a k - ir)PvU{a k ~ P w ,e 



w 



If k = 2s, and s — 1,3,..., then the probe state in the 
first term is Py/, while in the second term it is P&, and 
conversely if s is even. There remains (odd s) 

1 r> «. rr(l— *) r>(l — *;) T>(l—k) TT (l—k)^ 



+p*®u^- k) p { y- k) wp^- k) u {i - k) 



and so forth. 



C. Compensation 

Besides projections by P±, the fc-qubit system may 
experience some unitary transformations, as when 77 = 
— 1 for the qubit probe. One has 

u ; -k) u (i-k) J < _ u (i-k-i) u (v--k-i) J < p (k) 

2 TT (i-k-i) TT (l-k-t)i p (k) 



+\ft\ ur 'u'_ 

_ p(k) _|_ p(k) _ y 



using the recursion ( 3.8a ), and likewise for U-. It also 
follows that U± and P± commute. The unitary opera- 
tors U± are Hermitean, by recursive construction. Con- 
sequently 



{uf k) ) 



1 



(i-fe) 



Therefore, if the unitary transformations U± on the se- 
lected k qubits are not desirable, a clean projection can 
be obtained by repeating the measu rement, as described 
already for the Bell- aspect in ( 2.13 ). 

For the qubit probe case, there is an alternative way 
to compensate the unitary operations. For the k qubits, 
let 



uf z k) = ® • • • ® 



where 



This represents a multilateral +7r/2 rotation about the 
z-axis. Then it is straightforward to show that 



u i z 1 :: z k) u k ---u 1 

= U(a k ) ® g k P^- k) + U(0 k ) ® ig k P^ k) 
with 

9 = = ± (1 - i) 

For instance, one has 

U z P+=gP +1 U z P-=igP- 
Another way is to write 



(3.11) 



and use it for 9 = tt/2 as before, but with A 
the probe, i.e. 



P+ for 



U{9) 



e -iOP+ = 1 



-1)P+ 
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TABLE IV. GHZ-state propagation during two rounds of 
shuttle operation (overall phase omitted). The probe starts 
(and ends) in \E). Numbers refer to the listing in Table [n|. 



input 



first-round data 



output 



second-round data 



WESW 
WWSE 
EWSE 
EESW 
WENW 
WWNE 
EWNE 
EENW 



EWSE 
WWSE 
EWSE 
WWSE 
EWNE 
WWNE 
EWNE 
WWNE 



the qubits by U^, all in all (overbars indicate that the 
operators are built from a x instead of a z ) 

= {UlUsUy) {u}U 2 U y ) (ulUiUy) 

Finally the probe is measured in the y-direction, M y . 
Or, it is rotated by — 7r/2, say and measured with the 
M x equipmen t (t his option will not be included here). 
According to ( |3.6| ) the interactions produce 



£/(3tt/2) <8> U 



(123) p(123) 



+ U(w/2) ® U ( _>P_ 



(123) p(123) 



This gives r/ = U(—2tt) = 1. In the expression (3.6) one 
must then use U± = 1, together with 



U(a k ) = e~ ik 2 p + 



g k U(a k ) 



U(J3 k ) = e 



2 *■+ = ig*U{fi k ) 



Apart from the phase factors, the present global rotation 
compensates as effectively as repeating the measurement, 
although it requires a different set of operations from 
what is otherwise used, i.e. the U z . In some systems, 
this rotation could perhaps be drawn from an inherent 
'precession' of the qubits (not included here). 



D. Shuttle-measurement of the GHZ-states 



Here 



P (i23) = i (l(123)± ^3; ) 



(123h 



As the probe starts in Pe, say, when operating on the 
GHZ-states it will become either Pn or Pg , depending on 
the outcome of the e/o test, and hence must be measured 
in the y-direction. With rj = — 1, the unitary operators 
are 



U 



(123) 



ri 2 >Pi 3) 



P 



(3) 



r7 (123) =(J (2 )p (3 )+p (3) 



These determine which 3-qubit states are available after 
completion of the measurement. Use the following rep- 
resentation 



,(123) 



>(12) 



\(3) 



Entangled states with a structure such as the one 
shown in (3.2) for the eight GHZ-basis states can be mea- 
sured in a complete, non-demolition fashion by means 
of the partial measurement module just established. Of 
course, these modules can also be used separately in 
ways, which need not imply entanglement. In each stage, 
there are optio ns to repeat or to compensate, as described 
above in Sect. Ill C . For the following design, it will be 
regarded as acceptable that there are simple cycles be- 
tween the members of the GHZ-basis. 

The probe is initialized in the state Pe- It first in- 
teracts with qubits 1 and 2, by U% and {J 2 , and is then 
measured in the x-direction, M x . Secondly, it interacts 
with qubit 3 and again with qubit 2, via £/ 3 and L^- It 
is then measured again with M x . When the GHZ-basis 
states are given as input, this non-demolition processing 
leads to some permutations, as listed in Table III. They 
come from 



The first factor in the right-hand expression stands for 
any of the Bell-states, so 



00a 



00a 



^(12) = p(2) 



Pf = 1 

>(2) 



and the last factor is a a x eigenstate. The GHZ-basis in 
this representation is shown in Table ||. As usual, the 
even states 1 to 4 are operated upon by U+, while the 
probe turns 37r/2, and the odd states 5 to 8 by the 
probe turning ir/2. 

The resulting qubit output states and probe data are 
listed in Table III. It depends, of course, on the circum- 
stances, whether these output states are useful. In any 
event, there is a compensation procedure using only local 
rotations, U z , which will produce output within the origi- 
nal GHZ-aspect. Another such alternative output-aspect 
is entered if one omits the final . 

Repeating the last partial measurement corresponds to 
doing 

= {M x ul m u 1 u 2 u i u yyy ) (MJ : .l J .1 :i „„) 



c/ (12) = p(2)_p(2) =(j ( 2 ) ) 

C/| 32) = Pf - Pl 2) = , U m = Pi 2) + P?> = 1 

The pattern of probe data (2 bits) now allows to distin- 
guish four orthogonal subspaces, but there is not yet any 
certainty of entanglement. 

The e/o test starts by tilting each qubit spin by means 
oiUy. Then all qubits interact with the probe. Returning intermediate rotations by the unitarity of U yyy . At the 



According to Sect. IIIC, this will return the appropriate 
GHZ-states, i.e. the ones that are present after the first 



two stages in Table III 



This allows the following shuttle design. Cancel the 
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end, add Uy yy 



-iijyyy, which merely performs another 




Q, -h (N cn 



FIG. 2. Shuttle-operation of the complete, non-demolition 
measurement of the G HZ-a spect. The probe p interacts with 
qubits 1, 2, and 3 via (3.4) (wiggly lines) and is measured in 
the x- and y-directions at the standard measuring stations M x 
and My. The qubits are rotated by local U y . The GHZ-states 
pass unhindered, except for some permutations. 



permutation on the GHZ-basis (exchanging the even and 
odd states, i.e. (15) (26) (37) (48) apart from phase factors 
±1). This gives the shuttle procedure of Fig. || 

{M X UyyyU X V 2 V Z ) (MyU^U^yyy) 



The final output is shown in Table IV. One must of 



course do the third probe measurement, M y , in order 
to complete the projection of the qubit system into one- 
dimensional subspaces. The aspect at this stage is not 
GHZ. The data from the fourth probe measurement, M x , 
is redundant, but a useful check on the e/o test. The 
unitary operations preceding it produce the GHZ output- 
states listed. 

Having returned to the GHZ-aspect, the measurement 
sequence can now be repeated. Apart from phase-factors, 
this merely gives rise to the simple cycles 

• • • 1 1-> 5 i-> 1 • • • , • • -4 i-> 8 i-> 4- • • 

After two rounds, therefore, all the GHZ-basis states 
emerge in the original places, and the probe is in Pe, 
where it started. The GHZ-states also result when the 
input is an arbitrary mixed state, w, with statistical 
frequencies equal to the predicted probabilities for w. 
One can recognize the outcome of the complete, non- 
demolition measurement by the probe data sequence that 
has been recorded. For two rounds, as shown in Table IV, 
only eight 8-bit words out of a total of 256 possible ones 
correspond to an error- free measurement. These words 
are going to repeat as the shuttle continues running. 



IV. CONCLUSIONS 

The measurement procedure established in the present 
paper is different from the existing ones in at least one 
of the following respects: (a) It is a complete, non- 
demolition measurement process, under which all the en- 
tangled basis states are invariant, (b) The entire inter- 
action with the entangled system is taken care of by a 
single (qubit) probe, used repeatedly, (c) There is no in- 
teraction between the individual parts of the entangled 
system, (d) There is no interaction between the envi- 
ronment and the entangled system, only with the probe, 
(e) All interactions are two-qubit interactions between 
the probe and each of the constituents of the entangled 
system in turn, (f) This interaction is of the simplest 
conceivable form, and only one kind of probe-qubit in- 
teraction is required to perform the entire measurement, 
(g) The design is modular, so that it can perform partial 
measurements, and can be extended to handle entangled 
states of many qubits. (h) The same procedure can be 
adjusted to measure a wide range of different entangled 
states by local qubit operations, such as unilateral spin 
rotations. 



11 



Apart from this one needs external apparatus, which 
can perform measurements on the probe. It is not essen- 
tial how the apparatus is built, as long as it performs a 
standard measurement, say measures the probe spin in 
the x-direction. This equipment can also be localized, i.e. 
it does not require entangled states of the apparatus to 
be prepared and distributed beforehand. 

With the present procedure the measurement process 
is split up into exclusively local operations and interac- 
tions, which can take place sequentially as the probe vis- 
its each of the qubit components in turn. All these oper- 
ations are elementary, and in their totality they can be 
said to define the operational nature of the observable 
entanglement. 

It is possible to design the probe schedule in a shut- 
tle mode, which may be repeated indefinitely. In this 
way one can take advantage of the quantum Zeno effect 
to maintain the entangled states for a length of time, 
with an insignificant probability of dccohercncc. Such a 
device would be capable, in principle at least, to store 
entangled states until they are released at a predictable 
time, in a predictable state. The shuttle process is stable, 
and self-correcting, in the sense that any random errors 
are removed after a few steps by the probe measurement 
actions. 
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